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■ We investigate the stability of /(-R)-rotating (Kerr) black hole obtained from a 

limited form of f(R) gravity. In order to avoid the difficulty of handling fourth order 
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scalar equation leads to a massive scalaron equation. It turns out that the f(R)- 
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5-h ' black-hole bomb idea. 

(N 

t> . PACS numbers: 

o 

00 
CO 

q I Typeset Using ETgX 



X 



*e-mail address: ysmyung@inje.ac.kr 



1 



1 Introduction 



Modified gravity theories, f(R) gravities [H EJ [3] have much attentions as one of strong 
candidates for explaining the current accelerating universe [I]. f{R) gravities can be con- 
sidered as Einstein gravity with an additional scalar field. For example, the metric-/ (R) 
gravity is equivalent to the ojbd = Brans-Dicke (BD) theory with the potential [5]. It was 
shown that the equivalence principle test allows f(R) gravity models that are nearly in- 
distinguishable from the ACDM model in the background universe evolution [6] . However, 
this does not imply that there is no difference in the dynamics of perturbations [TJ. It is 
worth noting that the perturbation distinguishes between Einstein and f(R) gravities. 

In order for f(R) gravities to be acceptable, they obey certain minimal requirements for 
theoretical viability [21 E]- Three important requirements are included: (i) they possess the 
correct cosmological dynamics, (ii) they are free from instabilities, tachyon and ghosts [SJ E] , 
(iii) they attain the correct Newtonian and post-Newtonian limits. 

On the other hand, the Schwarzschild-de Sitter black hole was obtained for a constant 
curvature scalar from f(R) gravity [7J. A black hole solution was obtained from f(R) 
gravities by requiring the negative constant curvature scalar R = R [9]. For 1 + /'(_R ) > 0, 
this black hole is similar to the Schwarzschild-AdS (SAdS) black hole. In order to obtain 
the (constant curvature) black hole solution from f(R) gravity coupled to the matters of 
the Maxwell [9] and Yang-Mills fields [10], the trace of its stress-energy tensor should 
be zero. Interestingly, it was pointed out that the Kerr solution could be obtained from a 
limited form of f(R) gravity [TT]. Later on, it was argued that perturbed Kerr black hole 
can discriminate Einstein and f(R) gravities |12j . 

All astrophysical black holes belong to the rotating black hole. A black hole solution 
should be stable against the external perturbations because it stands as a physically realistic 
object [13] • Studies of stability of Kerr black hole are not as straightforward, because it 
is axially symmetric black hole and thus, the decoupling process seems to be complicated. 
The Kerr black hole has been proven to be stable against gravitational fields [HI HHJ HE] 
and massless scalar [TTJ. However, there exist unstable modes when considering a massive 
scalar due to the superradiance [TSj. Furthermore, it seems that the stability analysis of 
/(-R)-rotating black hole is a formidable task because f(R) gravity contains fourth order 
derivatives in the linearized equations [12] . 

In this work, we investigate the stability of f(R) -rotating (Kerr) black hole arisen from 
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a limited form (12. 6p of f(R) gravity [UJ. This work will be interesting because if f(R) 
gravity is considered as a viable theory, it is responsible to explain the present accelerating 
universe as well as the formation of rotating black holes. We transform the form of f(R) 
gravity into the scalar-tensor theory to avoid fourth order derivative terms by introducing 
two auxiliary scalars. Then, the linearized curvature scalar equation becomes a massive 
scalaron equation, indicating that all linearized equations are second order. Using the 
stability analysis of Kerr black hole in the massive Klein-Gordon equation, we show clearly 
that the /(i?)-rotating black hole is unstable against the scalaron perturbation. 



2 Perturbation of f(R) black holes 

We start with f(R) gravity without any matter fields whose action is given by 

S f = ^j d A xyf^f{R), (2.1) 
where k 2 = 8irG. The Einstein equation takes the form 

R,uf\R) - \g^f{R) + (^V 2 - V,V U ) f(R) = 0, (2.2) 



where ' denotes the differentiation with respect to its argument. It is well-known that (12. 2p 
has a solution with constant curvature scalar R = R. In this case, (12. 2p can be written as 

R,uf\R) - \g»uf(R) = o, (2.3) 

and thus, the trace of (12. 3p becomes 

R= 2 J^Sl = Ak f (2.4) 
f'{R) f 1 } 

with Af the cosmological constant due to the f(R) gravity. Substituting this expression 
into (12.31) . one obtains the Ricci tensor 

f(R) 

R »f = 2 f'(Rf flv = A/ ^' ( 2 ' 5 ' ) 

In order to find the Kerr black hole solution, we have to choose a non-pathologically func- 
tional form of f(R) as [TT] 

f(R)=a 1 R + a 2 R 2 + a 3 R 3 + --- (2.6) 
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which is surely a Talyor series around R = R = 0. We should mention that this form of 
f(R) gravity is not general, but a small subset of f(R) gravities. We check that /(0) = 
0, /'(0) = ai, /"(0) = 2a2 at R — R — 0, which provides either the Schwarzschild or Kerr 
black hole solution when choosing Af = 0. 

In this work, we use the Boyer-Lindquist coordinates to represent an axis-symmetric 
Kerr black hole solution with mass M and angular momentum J [TP] , 



/ 2Mr\ , 2 2Mrasin 2 # , ,, p 2 , , 
VP/ P A 



2 2 2Mra 2 sin 2 ^\ 
+p 2 d0 2 -( ( r + a H — ) sin 2 9 



(2.7) 



with 



A = r 2 + a 2 -2Mr, p 2 = r 2 + a 2 cos 2 9 , a = — . (2.8) 

In the limit of a — > 0, (12. 7p recovers the Schwarzschild black hole, while a — > 1 goes to the 
extremal Kerr black hole. The zeros of A, two horizons are located at 



r± = M ± VM 2 - a 2 (2.9) 
and the angular velocity at the event horizon takes the form 

n = —^—. (2.io) 

2Mr+ v ' 

Now we introduce the perturbation around the Kerr black hole to study stability of the 
black hole 

g^ = g f w + h fW . (2.11) 

Hereafter we denote the background quantities with the "overbar". Then, the linearized 
equation around /(-R)-rotating black hole becomes 



VVpKp + VV„/i w - V 2 V - V^V u h - g^(y a VPh a p - vV 

(^,V 2 - V M V„) (v Q V^ Q/3 - V 2 /i) = (2.12) 



3m 2 
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with the mass squared m 2 defined by 

= mry (2 - 13) 

Taking the trace of (I2.12p with g^, one has the fourth order equation for h^ u 

(v 2 - mf)8R(h) = -)• (V 2 - mj) (v a V' 3 ^ - V 2 /i) = 0. (2.14) 

At this stage, we note that it is not easy to make a further progress on the perturbation 
analysis because there exist fourth order derivatives. We mention that for the Einstein 
gravity with f(R) = R, /'(0) = 1 and f"(0) = 0. In this case, one finds the equation 
for linearized curvature scalar: SR(h) = 0, which means that SR(h) is not a physically 
propagating mode. Actually, this equation leads to one constraint as 

V a V p K p = V 2 h (2.15) 

which will also be recovered from the transverse gauge. 

Up to now, we did not fix any gauge. We would like to comment on the linearized 
equation by choosing the Lorentz gauge 



V„F = -V^h. (2.16) 
Under this gauge-fixing, the linearized equation (I2.12p takes a simple form 

V 2 /v + 2R ftpv<r h p<T + (g»u V 2 - V At V.) V 2 h = (2.17) 

with the trace-reversed perturbation h^ u = — and l/3m 2 = A [12]. This equation 

was importantly used to mention that perturbed Kerr black holes obtained from f(R) 
gravity can probe deviations from the Einstein gravity [12]. Even though equation (I2.17P is 
simpler than ( I2.12p . it is a non-trivial task to decouple odd and even perturbations around 
the Kerr black hole, arriving at two fourth order equations hopefully. Crucially, we do not 
know how to solve the fourth order equation (I2.17P around the /(-R)-rotating black hole. 
On the other hand, we may choose the transverse gauge which works well for studying the 
graviton propagations on the the Minkowski and AdS4 spacetime backgrounds [211 122] 



V„/i^ = Vh, (2.18) 
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which leads to (12.15P when operating V on both sides. Using the relation ( I2.15p . one 
immediately finds that the effect of f(R) gravity [2/3mj-term in (12.121) ] disappears because 
of 5R(h) = 0, leading to the Einstein gravity. Hence, the gauge- fixing for stability analysis 
of the /(-R)-rotating black hole differs from that for the graviton propagations on the AdS, 
dS, and Minkowski spacetimes [20]. We note that choosing an appropriate gauge-fixing 
cannot resolves the difficulty with f(R) gravities. It seems that the best way to resolve 
the difficulty confronting with the fourth order equation is to translate the fourth order 
equation into the second order equations by introducing auxiliary scalar fields. In other 
words, we must make a transformation from the limited form of f(R) gravity to the scalar- 
tensor theory (like Brans-Dicke theory) to perform the stability analysis of /(-R)-rotating 
black hole. 

3 Perturbation of the scalar-tensor theory 

In this section, we will develop the perturbation analysis around the /(-R)-rotating black 
holes (12. 7p in the different frame, the scalar-tensor theory. Introducing two auxiliary fields 
4> and A, one can rewrite the action (12 .ip as [231 [21] 

Sst = ^J d 4 x v ^[<f>(R-A) + f A (A)]. (3.1) 

Varying for the fields <fi and A leads to two equations 

R = A, <f> = f A (A). (3.2) 

Note that using (13. 2p . the action (13. ip recovers the original action (12.11) . On the other hand, 
the equation of motion for the metric tensor can be obtained by 

f' A (A)R, u - ^-g,u + (g»uV 2 - V,V u )f' A (A) = 0. (3.3) 

Considering a constant curvature scalar R = R = A together with = f' A (A) = const, 
Eq. (13.31) becomes 



f' A (A)R^ - l -g» u f A (A) = 0. (3.4) 



Taking the trace of (13 .4p leads to 
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Substituting this expression into (13. 4j) . one finds the Ricci tensor which determines the 
maximally symmetric Einstein spaces including Minkowski space 

At this stage, we introduce a specific form of /a(^4) inspired from (12. 6p [TT| 

f A (A) = ai A + a 2 A 2 + a 3 A 3 ■ ■ ■ , (3.7) 

where ai, a2, 03, • • ■ are arbitrary constants. Their mass dimensions are [ai] = 0, [a 2 ] = —2 
because [/^(A)] = 2, [A] = 2, and [0] = 2. In this work, we confine ourselves to the 
asymptotically flat spacetimes with = which accommodates the /(-R)-rotating black 
hole ( 12. 7p . In this case, we have to choose A = and thus, 

f A (0) = 0, f A (0) = ai > 0, f'i(0) = 2a 2 > 0. (3.8) 

Now we are in a position to study the perturbation around the Kerr black hole (12. 7p . In 
addition to ( 12. lip , from (13.21) . we have 

R + 5R(h) = A + 5 A, 4> + 5<fi = f' A (A) + f A (A)6A, (3.9) 

which leads to 

5R(h) -> 6 A, 5<f) -> f A (A)8A. (3.10) 

Thus, instead of 8R(h)-8(j), we use 5 A as a perturbed field in addition to h^ v . We mention 
that if one uses 5R(h)-5(p, one gets the same linearized equations. For this, see the appendix 
of Ref.[25]. 

Using (I3.10p . the linearized equations to ( 13. 3 j) reduces to 

SR^ih) = [^] V fi VJA+ ^JA, (3.11) 
where the scalaron mass squared is given by 

Taking the trace of (13. 1 1 j) and using (13.101) leads to the massive scalaron equation 

V 2 -m^W = 0. (3.13) 



Since the mass dimension of the linearized scalaron is two ([8 A] = 2), it would be better 
to write the canonically linearized equations by introducing a dimensionless scalaron 5 A 
defined by 

ft A 

M = ss- (3 ' 14) 



Finally, we arrive at two linearized equations 

V 2 -m^W = 0, (3.15) 



ftR^h) = V.VJA + V M, (3.16) 

which are our main equations to carry out the stability analysis of /(-R)-rotating black 
hole. It is a nontrivial task to perform the stability analysis of /(-R)-rotating black hole 
based on the metric-perturbing equation (13 . 16[) since the decoupling process seems to be 
complicated. However, we expect that the right hand side of scalaron coupling terms in 
( I3.16P will not change the stability of f(R) -rotating black hole (23 |2B]. In this work, we 
focus on the scalaron equation (13.15P because it shows the feature of f(R) gravities. 



4 Stability analysis of /(i^)-rotating black hole 



Considering the axis-symmetric background (12 .7p . it is convenient to separate the scalaron 
into [27] 



5A(t,r,6, 



im ^S] n {e)R{r) 



(4.1) 



where S™^) are spheroidal angular functions, and the azimuthal number m takes on (posi- 
tive or negative) integer values. Also, it is enough to consider positive cj's in (14. ip . Plugging 
this into the scalaron (Klein-Gordon) equation (I3.15p . we get the angular and radial wave 
equations for S™{6) and R{r) as 



sin^ 



deismOdeS? 



a 2 (oj 2 — m 2 A ) cos 2 9 



m 



sm z 6 



A 



lm 



ST- 



A<9 r (Ad r R) + [u 2 (r 2 + a 2 ) 2 - AMamur + a 2 m 2 
-A(a 2 u 2 + m 2 A r 2 + A lm )] R = 



(4.2) 
(4.3) 
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where Ai m is the separation constant that allows the split of the wave equation [28j 129] 

oo 

A lm = l(l+l) + J2 c k a 2k (m\ - cu 2 ) k . (4.4) 

k=l 

The radial Teukolsky equation takes the Schrodinger form 

d 2 R 



, - + V(r, u)R = u 2 R, R = Vr 2 + a 2 R, (4.5) 
dr* 2 

where the tortoise r* coordinate is defined by dr* = r2j ^f dr and the w-dependent potential 
V(r,cu) is given by 

Am 2 , AMramu — a 2 m 2 + A\Ai m + (u 2 — m 2 A )a 2 ] 
V(r,u) = 2 , 2 + 



r 2 + a 2 (r 2 + a 2 ) 2 

A(3r 2 - 4Mr + a 2 ) 3A 2 r 2 
(r 2 + a 2 ) 3 (r 2 + a 2 ) z 



^ ~ — ' ^ ; _ — (46) 



Its asymptotic forms are given by 

V -> u 2 - m 2 A , r* -> oo(r -> oo), (4.7) 

V-^(w-mfi) 2 , r* -> -oo(r -> r+). (4.8) 
Near the horizon at r = r + and spatial infinity at r = oo, the scalaron takes the form 

R = T e -K<»-™£l)r* ; r * _ OQ (4 9) 

£ = e-'^^' + ^V" 2 -**, r*^oo (4.10) 

with the T(TZ) the transmission (reflection) amplitudes. Here we must assume that u, for 
the bound case of cu < uia, lies on the physical sheet of < arg(w 2 — m 2 ^) 1 / 2 < ir. Requiring 
ingoing waves upon a rotating black hole whose angular velocity Q is given by (I2.10p . one 
must impose an ingoing (negative) group velocity v gr for the wave packet. Then, we choose 
an ingoing mode near the even horizon as 

[e-^Jin = Te- iwt e- i{uj - mn)r * . (4.11) 

Since V(r,u) is real, the Wrongskian W(R,R*) of the complex conjugate solutions of R 
and R* satisfies [13] 

ijL W (R,R*) = (4.12) 
dr* 
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which implies that 



n\ 2 = 1 + 



mVt 



i in 2 . 



(4.13) 



Here, if the frequency oo of the incident wave satisfies the inequality 



oo < mVt , 



(4.14) 



one has \1Z\ > 1, implying that the reflected (scattered) wave is being amplified j27J HSj- 
Thus, in this superradiance regime, waves appear as outgoing to an observer at infinity, and 
energy radiated away to infinity actually exceed the energy present in the initial perturba- 
tion. Feeding back the amplified scattered waves, one may gradually extract the rotational 
energy of the /(i?)-rotating black hole. Press and Teukolsky have suggested to use this 
mechanism to the black-hole bomb [13]. If one surrounds the black hole by a reflecting 
mirror, the wave will bounce back and forth between black hole and mirror, amplifying 
itself each time. To this end, nature may provide its own mirror [281 130] . 

It is well known that if one considers a massive scalar field with mass /i scattered off a 
rotating black hole, then for u < fi, the superradiance has unstable modes because the mass 
term effectively works as a reflecting mirror [311 Q2J E21 H31 EH HE] • As was shown in Fig. 15 
of Ref.|13]. a shape (ergoregion-~-mirror) of potential V(r, u) has the local maximum as well 
as the local minimum far from the black hole which generates a secondary reflection of the 
wave the reflected from the potential barrier. This implies that the meta-stable bound states 
of Mm a ~ r + /X c with A c the Compton wavelength of the massive particle could be formed 
in the valley of the local minimum. The secondary reflected wave will be reflected again at 
the far region. Since each scattering off the barrier in the superradiant region increases the 
amplitude of the wave, the process of reflections will continue with the increased energies of 
waves, leading to an instability. In our case, the mass tua of scalaron works as a reflecting 
mirror. Considering a massive scalaron with mass tua scattered off /(-R)-rotating black 
hole, then for u < tra, the mass term effectively works as a reflecting mirror. Similarly, 
the Kerr black hole was shown to be unstable when choosing f(R) = R + hR 2 [3H] which 
is included as a limited form ( 12. 6 p when setting all a n = for n > 3, a\ — 1, and ci2 = h. 

The instability time scale r associated with the dynamics of a massive scalar was re- 
stricted to two limiting cases: Mvtla 1 provides r ~ 24(a/M) _1 (Mm^)~ 9 (GM/c 3 ) [IS] , 
while Mm A > 1 indicates r ~ 10 7 e L84MmA (GM/c 3 ) [32]. Recently, the authors [M] have 
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shown that for Mm^ < 0.42 and a = 0.99, the maximal instability growth rate is given by 
t- 1 ~ 1.5 x 10~ 7 (GM/c 3 )' 1 . 

5 Discussions 

It was suggested that perturbed Kerr black holes discriminates Einstein and f(R) gravi- 
ties [12] . To this end, we have made a significant progress if the scalaron approach (scalar- 
tensor theory) represents truly the f(R) gravity including the fourth order differential 
equation. 

It turned out that the f(R) black hole is stable against external perturbations if the 
scalaron does not have a tachyonic mass {m\ = /a(0)/3/^(0) > 0, /^(0) > 0) [25]. This is 
consistent with other perturbation analysis: the Dolgov-Kawasaki instability with f"(R) < 
in cosmological perturbations [37]. For the f(R)-AdS black hole, its stability is guaranteed 
if the scalaron mass squared rh\ satisfies the Breitenlohner-Freedman bound (rh\ > m^ ¥ = 
—9/4£ 2 ) [26]. However, for /(i?)-rotating black hole, it is unstable even for the positive 
scalaron mass squared m\ > because of its superradiant instability (black hole bomb). 
We would like to stress that this is a meaningful result for the rotating black hole obtained 
from a limited form (12. 6p of f(R) gravity. 

Given the results for the e-folding time r = 10 7 e 1MMmA (GM/c 3 ) for Mm^ ^> 1, when 
might the instability be significant rTSJ E5]? One assumes that if the scalaron is a pion 
around a solar-mass black hole, MrriA ~ 10 18 with M ~ 10 30 kg. Then, one concludes that 
the instability growth rate r _1 is not significant for astrophysical black holes, unless there 
exists an unknown particles with a tiny but nonzero rest mass. That is, the instability 
growth rate r _1 is always small for the standard model particles, in compared to the decay 
rate of particle or the Hawking evaporation rate of the black hole. However, the instability 
may be important for a pion around a small primordial black hole (M < 10 12 kg). 

Finally, f(R) theories are supposed to be viable only when the matter is present outside 
the star/black hole so that the Chameleon mechanism could take place. However, in order 
to set matter (even a small amount of matter) outside the star/black hole, it will change 
the results of analysis for the f(R) theories, as the scalaron 5 A becomes massive. 

In conclusion, we have shown that the /(-R)-rotating black hole arisen from a limited 
form (12. 6p of f(R) gravity is unstable because of the superradiant instability known as the 
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black-hole bomb idea. 
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